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We study how the electron hopping reduces the Mott-Hubbard band gap in the limit of a large 
Coulomb interaction U and as a function of the orbital degeneracy N. The results support the 
conclusion that the hopping contribution grows as roughly i/NW, where W is the one-particle band 
width, but in certain models a crossover to a ~ NW behavior is found for sufficiently large N. 
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I. INTRODUCTION 

The Mott-Hubbard metal-insulator transition has aU 
tracted much interest since it was proposed by Mott,EJ 
being an example of how strong correlation drastically 
changes the properties of the system. For simplicity, the 
Mott-Hubbard transition is-,usually studied for models, 
without orbital degeneracyH e.g., the Hubbard model.cl 
For most real systems, however, the orbitals primarily in- 
volved in the transition are degenerate. We recently ar- 
gued that in a Hubbard model with an orbital degeneracy 
N, the ratio of U/W where the Mott-Hubbard transition 
takes place is increased by roughly a factor V^V,Q where 
U is the on-site Coulomb interaction and W is the one- 
electron band width. 

This conclusion was based on a simple and suggestive 
but nonrigorous argument in the large [/-limit, which 
said that the band gap E g behaves aaj 



E a ~U- VNW. 



(1) 



This argument was supported by exact diagonalization 
calculations for small clusters with N = 1, 2 and 3. The 
argument was then extrapolated to intermediate values 
of U, and the extrapolation was supported by quantum 
lattice Monte Carlo calculations for a model of A3C60 
(A=K, Rb) with the N = 3. Similar results have recently 
also been-,obtained in Monte Carlo calculations by Han 
and Cox,cl while Lrj3 found a (N + 1)W behavior using 
the Gutzwillcr Ansatz and Gutzwiller approximation. 

Here we want to investigate the behavior in the large 
U limit further. We have in particular studied a model 
with hopping only between orbitals with the same orbital 
quantum number, referred to as an intraband hopping 
model. This model shows an interesting cross-over as the 
degeneracy N becomes larger than the number of nearest 
neighbor K. Thus the hopping contribution tends to be 
- VNW for N « K but - NW for N » K. When 
hopping between orbitals with different quantum num- 
bers is included, interband hopping models, the results 
become more model dependent. The results obtained 
here support a behavior ~ >/NW or a somewhat more 
rapid growth with N. 



In Sec. II we present some general arguments. In Sec. 
Ill we use exact diagonalization to study small clusters 
with intraband hopping and in Sec. IV this is extended 
to include interband hopping. The results are discussed 
in Sec. V. 



II. GENERAL CONSIDERATIONS 

We consider a system with M sites and the orbital de- 
generacy N at half- filling, i.e., with NM electrons. The 
band gap is then given by 

E g = E(NM + 1) + E(NM - 1) - 2E(NM), (2) 

where E(L) is the ground-state energy for L electrons. 
In the large (7-limit, the ./VL-electron ground-state has 
exactly N electrons per site. The hopping of an electron 
to a neighboring site would cost the Coulomb energy U 
and is strongly suppressed in this limit. Therefore 



1 ft 2 
E(NM) = -N(N - l)MU + O ( — 



(3) 



where t is a hopping integral. For the states with an extra 
electron or hole the situation is more complicated, since 
the extra electron or hole can hop without any extra cost 
in Coulomb energy. We therefore focus on these states. 
We first form an antiferromagnetic state |anti) with N 
electrons of a given spin on each site and with each site 
having as many neighboring sites of the opposite spin as 
possible. To be specific, the cental site 1 has spin up. We 
then form a state with NM + 1 electrons 



K) = ^ni |anti), 



(4) 



where ip\ ma creates an electron on site i in the orbital m 
with spin a. Thus site 1 in the state \v\} has an addi- 
tional electron. A spin up electron can hop from this site 
to a neighboring spin down site, forming a state 
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V4uXMmT^ lm ~rl anti ); 



(5) 



where the sum i is over the L (< K) neighboring sites 
with opposite spin to site 1. The spin down electron on 
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site 1 can also hop to neighboring sites with spin up, 
but we neglect this for the moment. The matrix element 
between the states in Eqs. (0-pJ) can then be calculated. 



(V2\H\vi) 



(6) 



where we have considered the case when there is a hop- 
ping matrix element only between orbitals of the same 
m quantum number on the neighboring sites. The a&, 
gument can also be repeated for more general hopping.cl 
The important aspect of the result @ is that it has a 
factor y/N relative to the one-particle case.oQ This is a 
consequence of the fact that any of the N spin up elec- 
trons can hop to the neighboring site, as illustrated in 
Fig. [j]. For a bipartite system with N = 1, the ferromag- 
netic arrangement in Fig. [l]a is favorable according to 
Nagaoka's theoremB For N > 1 the situation in Fig. [l]b 
is in general more favorable, since the hopping can then 
take place in N channels. 



a) 



AAA 



AAA 



b) 



AAA 









\ 


1 


T 



FIG. 1. Illustration how an extra electron can hop against 
a half-filled background, a) ferromagnetically aligned neigh- 
bor: There is only one hopping-channel. b) antiferromag- 
netically aligned neighbor: There are N different ways for 
hopping. 

For one single electron in the system, the hopping en- 
ergy is of the order of W/2, where W is the one-electron 
band width. The extra factor y/~N in Eq. (|^) suggests 
that the hopping energy is correspondingly larger for the 
NM + 1 electron case and that the total energy is 



(7) 



E(NM + 1) » E(NM) +NU + VNW/2 



Analogous arguments then suggest a similar result for 
E(NM - 1) and we obtain 



E a = U 



NW. 



(8) 



As pointed out in Ref. [|, these arguments are not rig- 
orous. As the extra occupancy moves through the system 
it reduces the spins along its path. There is therefore no 
one to one correspondence to the one-particle case. We 
therefore need to consider the arguments in more detail. 

To be specific, we consider a degenerate multiband 
Hubbard model 



^im.jra V 'im CT V 'jm 'a 



(9) 



with hopping only between orbitals on nearest neighbor 
sites. In part of the discussion we furthermore consider 
the intraband hopping model, where t im Jm ' = t5 mrn > , 
i.e., only hopping between orbitals with the same quan- 
tum number m. 

To calculate the energies of the states with ran extra 
hole or electron, we use the Lanczos methods In this 
method a set of basis states \v n ) are generated according 
to the prescription 

b n+1 \v n+1 ) = H\v n ) - a n \v n ) - b n \v n -i) (10) 
a n = (v n \H\v n ) 

where b n+ i is determined by the requirement that |«n+i) 
is normalized. In this way the Hamiltonian is trans- 
formed into a tridiagonal form, where the diagonal el- 
ements are given by a n and the nondiagonal elements by 
b„ (H rm = a n and if„,„+i = b n+ i). The starting state 
is defined in Eq. (Q) and |u ) = 0. 

This approach gives the energy of the lowest eigen- 
state which is nonorthogonal to the state \v\). If our 
choice of \v\) is orthogonal to the ground-state of the 
(NM ± 1)— system, the energies of these ground-states 
are overestimated. We then obtain an upper limit to E g , 
since we know the energy of the half-filled state in the 
large i7-limit. This remains true if the iterative steps 
implied by Eq. ( |l0|) are interrupted before convergence, 
since the lowest eigenvalue obtained from the tridiago- 
nal matrix is decreased as more (nonzero) elements are 
added to the matrix. We furthermore observe that if 
\vi) is chosen poorly in the sense that it is very different 
from the ground-state, the ground-state energy is still 
obtained after sufficiently many steps, as long as \vi) is 
not orthogonal to the ground-state. 

In the following we often use a classical Neel state to 
construct the states |anti) and \vi). These states are con- 
vienent to handle, but not necessarily the best starting 
point, due to the tendency of the system to form a singlet 
(or a doublet with an extra electron or hole for N > 1). 
Due to the arguments above this will at worst mean that 
we underestimate the hopping reduction of the gap. 

In the case of a half-filled system with an additional 
hole or electron, the Coulomb energy is minimized if there 
are exactly N electrons per site, except for one site with 
the extra electron or hole. Since we consider the large- 
U limit we allow only such states in the following. The 
Coulomb energy is then a constant which we choose as 
energy-zero. We thus arrive at a t-J model with J = 0. 
Let T lo (NM ± 1) and T U (NM ± 1) denote the lowest 
and highest eigenenergy of the Hamiltonian. In the large 
{/-limit 
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T lo (NM ± 1) = -T hi (NM T !)■ (11) 

It is therefore sufficient to study, e.g., the (NM + 1)- 
electron case, which gives 

E g = U + T lo (NM + 1 ) - T hi (NM + I). (12) 
III. INTRABAND HOPPING 



From the results above it is immediately clear that for 
sufficiently large values of N, the reduction of the band 
gap due to hopping must at least be of the order Nt even 
for a larger system. This can bee seen by simply, arbi- 
trarily, picking two neighboring atoms in the system of 
interest and then constructing a state like above. This 
already gives a reduction of the band gap ~ Nt ~ NW, 
and any improvement of this simple construction can only 
make the reduction larger. 



We have studied a few small clusters of atoms, where 
the system is small enough to allow an exact diagonal- 
ization. We first consider intraband hopping models. In 
these models the one-particle band width of each channel 
m is the same and independent of N. 



A. Diatomic molecule 

We first study the case of just two sites, orbital degen- 
eracy N and 2N + 1 electrons. We form one state \a), 
with N + 1 electrons on site 1. This state is a linear com- 
bination of states, where the first state has N spin up 
electron and one spin down electron with, say, m — 1 on 
site 1 and N spin down electrons on site 2. To generate 
further states in \a), we let one of the spin up electrons 
hop to site 2 and one spin down electron on site 2 hop 
to site 1. This gives N(N — 1) states which are added 
to state | a). In this way we form all possible states with 
N + l electrons on site 1 and include them in \a). In this 
process we also generate all possible states with N + 1 
electrons on site 2, which are included in a state |6). For 
symmetry reasons, state \a) and \b) include equally many 
terms L and therefore have the same normalization factor 
1/y/L. 

In calculating the matrix element (b\H\a), we notice 
that in each (m,a) channel except (m = 1,|), there 
is only one electron. Each one of the L states in \a) 
then connects to exactly N states in \b), since N of the 
electrons (spin up or down) on site 1 can hop to site 2, 
and only the electron in the orbital (m = 1,1) cannot 
hop. There are then NL contributions t to the matrix 
element. After taking the normalization factor into ac- 
count, we then find that the matrix element is Nt. The 
corresponding result for the gap in the large U limit is 
then 



Eg = U - 2Nt = U - NW. 



(13) 



The reduction of the band gap due to hopping is thus 
proportional to N and not y/N as the arguments in the 
previous section suggested. If we start the Lanczos pro- 
cedure from the antiferromagnetic state discussed before 
(Eq. (§)), b 2 = VNt. For n ~ N, however, b n ~ Nt. 
The ground-state of the (27V + l)-electron system is a 
doublet. It is then not surprising that a large number 
of Lanczos steps is needed to generate this state from 
antiferromagnetic state \v\) in Eq. (0). 



B. "Bethe" lattice 

To obtain a system where an atom has several neigh- 
bors, we have studied a simple version of a finite Bethe 
lattice, where atom 1 connects to K neighboring atoms, 
while these atoms only connect back to the first atom 
and has no other neighbors. We have solved this prob- 
lem using exact diagonalization for a number of values of 
N and K and find that the results are described by 



Eg = \U-y/N + W 

1 U - NW 



1)/KW ifiT>l; 

if K = 1. 



(14) 



This result shows an interesting cross over behavior. For 
N << K, the reduction of the band gap is proportional 
to \ZWW, as suggested by the argument in Sec. II. 
On the other hand, for N >> K the reduction is pro- 
portional to NW, as suggested by the previous section 
(IIIA). The limit N < K should apply to most cases 
of practical interest, but the result for N >> K is also 
interesting. 



a) o- 
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FIG. 2. Schematic picture of small clusters studied by 
exact diagonalization. Each circle indicates a site and the 
lines show the hopping between these sites. 
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C. Other small clusters 

The results for the Bethe cluster are instructive, but 
because of the simplicity of the model it is not clear to 
what extent they apply to more realistic models. We have 
therefore studied a few models which are small enough 
to allow exact diagonalization. In Fig. || we show these 
models schematically. We define a quantity 



C(N) = 



U-E a (N) 



1 



lim 



U^ooU-E g (N = 1) v/jV' 



(15) 



where E g (N) is the band gap for the orbital degener- 
acy N. C(N) = 1 implies that hopping reduction of the 
band gap is proportional to yN. The results are shown 
in Fig. g. The figure supports the results found in Sec. Ill 
for the "Bethe" lattice, in the sense that C(N) is closer 
to a constant the larger K is, but that C(N) grows when 
JV~ K. 



U 




a (K = 2} 
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FIG. 3. The ratio C(N) between U — E g for a given value 
of the orbital degeneracy N and the same quantity for N = 1 
divided by y/N. Thus if the reduction of the band gap were 
proportional to i/N, C(N) would be identically unity. The 
figure illustrates how this tends to be true when the number 
K of nearest neighbors is much larger than N. The labels re- 
fer to Fig. ^| We have also included the results for a diatomic 
molecule (K = 1). 



We notice that the band width is reduced by the frustra- 
tion, since it is not possible to find a fully anti-bonding 
state. 

It is instructive to obtain this result from the Lanczos 
method. The nonzero coefficients a n and b n are shown in 
Table |. The process stops already after two step, since 



H\V2) 



0,2 m) 



b 2 \vi)=0. 



(17) 



We note that in the present case an electron on a given 
site can hop to four other sites. The effect of one such 
hop is cancelled by the term b n \v n —i). The remaining 
hops may contribute to 6 n +i or to a n . In the present 
case the contribution goes to b n +i in the first step but 
to a n in the second step. For the many-body case the 
coefficients a n are typically smaller and a larger part of 
the hopping contribution instead goes to b n +i. 

For N = 1, 2 and 3 the lowest possible (NM ± 1)- 
energy is reached after just a few (4-8) Lanczos steps, 
illustrating that the choice in Eq. (Q) of \vi) leads to the 
ground-states. The situation is different for N — oo. We 
have calculated the first 40 coefficients b n +i, of which 
the first few are shown in Table |. At least up to n = 40 
bn+i/VN keeps increasing with n. This suggests that the 
reduction of the band gap due to hopping grows faster 
that y/~N, in agreement with the arguments in Sec. III.A- 
B that a contribution ~ N may be expected for A >> K . 

TABLE I. The coefficients a n and b n for a tetrahedron. 
Both the one-particle and (NM — 1) many-particle problems 
are considered, where in the latter case the orbital degenera- 
cies N = 1, 2, 3 and oo are considered. The coefficients are 
given in units of ts/~N , where t is the hopping matrix element. 



n One-particle 



Many-particle 



N = 1 N = 2 N = 3 iV = oo 

0. n &n + l «n &n + l dn &n+l 1 n 6 n + l Ctrl ?>n+l 



0.00 
2.00 



1.73 
0.00 



0.00 
0.67 
0.19 
0.25 
-0.11 



1.73 
1.25 
1.51 
1.07 
0.00 



0.00 
0.57 
0.44 
0.31 
-0.17 
1.15 
-0.17 



1.58 
1.30 
1.48 
1.21 
0.98 
0.97 
0.00 



0.00 
0.49 
0.45 
0.24 
-0.00 
0.55 
-0.20 
0.94 
-0.16 



1.53 
1.32 
1.49 
1.34 
1.37 
1.21 
0.92 
0.91 
0.00 



0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



1.41 
1.41 
1.58 
1.76 
1.94 
1.96 
2.24 
2.17 
2.46 



D. Tetrahedron 

In this section we study the tetrahedron ( c) in Fig. 
in somewhat more detail. This system has four sites 
with equal hopping matrix elements t < between all the 
sites. In the one-particle case, there is one nondegenerate 
eigenvalue 3t and a three-fold degenerate eigenvalue -t. 
The band width is then 

W = At. (16) 



To discuss the large AT- limit, we first discuss the gen- 
eral situation when each site is surrounded by L{< K) 
other sites with the opposite spin. To leading order in N 
there are then L ways of hopping to a neighboring site 
increasing the number of flipped spins. In addition a spin 
flipped in a previous hop can be restored. This is immedi- 
ately seen by applying the Hamiltonian to \i)2) in Eq. (|5|) 
which can take us back to \v\). This latter contribution 
is cancelled by the term — b n \v n -\) in Eq. ( |To[ ) in the 
first few steps. Furthermore a n /\^N = in the N — > oo 
limit, for any finite n. We then find that 6„+i = V LNt 
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for the first few terms. This is illustrated in Table Q, 
where L — 2. 



If all the elements b n+ i would stay constant in the large 
AT- limit, b n+ i — a*J~Nt (a > 0), it is easy to show that 
the lowest state has the hopping energy — 2a^J~Nt and 
the band gap is E g = U — Au^fNt for large U. In reality, 
however, after a few steps the situation becomes more 
complicated than discussed above. 

After three steps, terms proportional to 



V4m T ^ lm T^ m - i V'2 m '|V'2 m " T , /'i m "Tl anti )' ( 18 ) 

mm m 




enter the calculation. Since the states with m < m are 
equal to the states with m > m , the norm of the state 
in Eq. (|f) is 4N 2 (N _ l)/2 « 27V 3 . As a result b 3 in 
Table | is y / 2.5Nt instead of y/2Nt. 

After four steps certain states can be reached in two 
different ways. This is illustrated schematically in Fig. ||. 
The amplitudes for these states then have to be added. 
The contribution of such a state to & 2 +1 is then not 
~ (l 2 + l 2 ) as assumed above but ~ 2 2 . This increases 
the value of & ra +i . In addition, the restoration of a flipped 
spin can be done in different ways, leading to new states 
not contained in b n \v n — i) in Eq. JlO|). This is illustrated 
in Fig. n, and it also increases b 2 +1 . 



(d) 




FIG. 4. The arrows illustrate how the double extra occu- 
pancy hops from site 1 to 2 and back, followed by a hop to 
3 and back. In each hop a spin is reduced. The same final 
state can be obtained by hops 1— » 3 — * 1 — ► 2 — ► 1, and the 
corresponding amplitudes are added. 



FIG. 5. a) Shows one contribution to \v n ). b) and c) 
show contributions to H\v n ). b) Shows how a contribution 
in |« n -i) can be recovered by allowing the extra occupancy 
to hop back from site 2 to site 3 and thereby increasing the 
spins of sites 2 and 3. c) Shows how the extra occupancy 
instead can hop to site 1, increasing the spins of site 1 and 2 
and creating a new state in |v„+i)- 

E. Bipartite and non-bipartite systems 

A bipartite system is a system where it is possible to 
partition the system in two sublattices with hopping only 
between the sublattices. Nagaoka's theoremH states that 
for a such a system with the orbital degeneracy N = 1 
and U = oo, the system is ferromagnetic when it has one 
electron or one hole relative to half-filling. The electron 
or the hole then have the same hopping possibilities as 
in the one-electron case and it follows that in the large 
{/-limit Eg — U — W, with the prefactor exactly equal 
to one. For N > 1 or for many of the cases with N = 1 
but a non-bipartite lattice, where the Nagaoka's theorem 
has not been proven, the reduction of E g due to hopping 
is larger than W. This illustrates that Nagaoka's theo- 
rem is not valid in these cases. Cases with N = 1 but a 
non-bipartite lattice are studied below. 

For a bipartite lattice where all nonzero hopping in- 
tegrals are equal, we can always construct states which 
are fully bonding or antibonding, by letting the coeffi- 
cients of the orbitals on the two sublattices having the 
same or the opposite sign. This is in general not possible 
for a non-bipartite lattice. For instance, for the tetrahe- 
dron we found a bonding state with the energy 3t but 
the antibonding states have only the energy —t. In the 
Lanczos formalism, this result comes about because an 
electron can hop in a closed loop using an odd number of 
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hops. This leads to nonzero coefficients a n , which makes 
the spectrum nonsymmetric around e — 0. This also re- 
duces the coefficients b n+ i as discussed in Sec. III.D. In 
the many-body case, the hopping of an electron along a 
closed loop usually does not bring us back to the original 
state, since spins have been flipped along the path of the 
electron. This tends to reduce the values of a n and to 
increase b n , as is illustrated in Table |[ Even for N = 1 , 
hopping may therefore reduce the band gap more than 
one would expect from the one-particle band width. This 
is illustrated in Table ||. In the cases we have studied 
the effect is of the order 15-30 %. 

TABLE II. The reduction U - E g due to hopping for 
N — 1 relative to the one-particle band width W for some 
non-bipartite systems. The notations of the systems refer 
to Fig. ^ and "triang" refers to a system with three atoms 
coupling to each others. 



System 


[U-E g (N = l)]/W 


c 


1.25 


d 


1.15 


c 


1.19 


triang 


1.33 



IV. INTERBAND HOPPING 

In the previous sections, we have for simplicity as- 
sumed that there is only hopping between orbitals on 
different sites with the same m quantum number. Here 
we consider the more general case of hopping also be- 
tween different m quantum numbers. We first consider 
the step in the Lanczos procedure. In the one-particle 
case, this gives 



(& 2 m) 



E k, 



(19) 



where the electron is located in orbital m on site 1 in the 
starting state. In the many-body case, we start from the 
state in Eq. (0) and obtain 



(M 2 = E 



(20) 



We compare this with the average of the coefficients 
(& 2 m) ) 2 in the one-particle case 



(21) 



i.e., 62 is typically a factor y/N larger in the many-body 
case than in the one-particle case. This is the same fac- 
tor as was found when the hopping between different m- 
quantum numbers was neglected. The inter-orbital hop- 
ping increases 62 by a factor y/~N in both the one-particle 



and many-body cases, but leaves their ratio (= y/~N) un- 
changed. 

TABLE III. The band gap for a six-atom cluster ( d) in 
Fig. ^). In all cases the absolute value of the hopping matrix 
elements is t. In the "sp" -model we have chosen the signs as 
if the two orbitals had s- (N — 1) and p- character (N = 2). 
In the "rand" model the signs were chosen randomly. For 
this model K = 4 and for the "sp" -model K is effectively 
somewhat smaller. 



Syst 


d(N) = [U- 


E g }/[W(N)VN] 


d(2)/d(l) 




N = 1 


N = 2 




"sp" 


1.15 


1.31 


1.14 


"rand" 


1.23 


1.10 


0.89 



In Eqs. (H© the sign of t im - m ' does not matter. 
When we continue to higher coefficients b n the sign be- 
comes important and we have to specify the model in 
more detail. We can, for instance, specify the orbitals 
involved and the geometrical structure of the lattice. In 
Table HI we show results for the lattice d) in Fig. |^, in- 
cluding a s-orbital and a p-orbital with the lobe along 
the direction of the three atoms. This defines the signs 
of the hopping integrals. The magnitudes are assumed 
to be identical. Alternatively, we can pick the sign ran- 
domly. The latter procedure is appropriate when we want 
to vary N, since it is not clear how to define an appropri- 
ate set of orbitals for an arbitrary value of N. Table III 
shows the ratio 



d(N) 



U-E g (N) 



W(N)VN 



(22) 



where W(N) is the one-particle band width for the or- 
bital degeneracy N. In contrast to the intraband model, 
the one-particle band W(N) width now depends on N, 
We also show the ratio d(2)/<i(l), which would be unity 
if the reduction of band gap due to hopping were propor- 
tional to VNW(N) . This is roughly the case in Table [n| 
but the results are too limited to allow more definite con- 
clusions. 

For the diatomic molecule, discussed in Sec. III. A, and 
the "Bethe" lattice, discussed in Sec. III.B, we can go to 
larger values of N. The results are shown in Fig. ||. 
These results support a reduction of the band gap due 
to hopping which is at least proportional to y/NW(N). 
Although the results are still too limited to draw definite 
conclusions, they do not support a behavior proportional 
to NW(N). The NW(N) behaviour found for the di- 
atomic molecule with only intraband hopping, is there- 
fore probably an artifact of the simple model used in that 
case. 
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g, 

T3 



1.45 

1.4 
1.35 

1.3 
1.25 

1.2 
1.15 

1.1 
1.05 
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FIG. 6. The ratio d(N) (Eq. (|22j) between the reduction 
of the band gap due to hopping and y~NW(N). The fig- 
ure shows results for the diatomic molecule (K — 1) and the 
"Bethe" lattice with K = 2 and if = 3. 



V. DISCUSSION 

Above we have studied two types of models, namely 
models with only intraband hopping (between orbitals 
with the same quantum number m) and models also with 
interband hopping. We have calculated the reduction of 
the band gap due to hopping in the large [/-limit for 
different orbital degeneracies N and for lattices with dif- 
ferent number K of nearest neighbors. In the intraband 
models we find that the hopping contribution grows at 
least as y/NW, where W is the one-particle band width. 
For N > K there is an interesting crossover, and the 
hopping contribution grows as NW . This limit should 
not, however, apply to most systems of interest. In the 
interband models the hopping contribution seems to be 
closer to a \/~NW(N) behavior than a NW(N) behavior 
even when K is small. Our present results therefore sup- 
port a degeneracy dependence of roughly ^/NW(N) or a 
somewhat stronger N dependence for most systems. 

It has been suggested that in a strongly correlated sys- 
tem the hopping should be strongly suppressed and it 
should not reduce the band gap much. For instance, in 
models of the High T c compounds, e.g, the t-J model, it 
is found that the string of flipped spins in the trace of a 
moving particle tends to strongly reduce the dispersion of 
this particle and that in addition the weight of the quasi- 
particle is strongly reduced. This might suggest that hop- 
ping is very inefficient in this case. In the present prob- 
lem we are, however, not interested in the quasi-particle 
dispersion but in total energies (Eq. (0)). The total en- 
ergy can be obtained by integrating the Green's function 
over frequencies and momenta£3 Since the weight of the 
quasi-particle is strongly reduced in correlated systems, 
this means that the spectral function must have substan- 
tial weight also for other frequencies. To argue in terms 
of the quasi-particle alone neglects important contribu- 



tions to the total energy. It is then more convenient to 
follow the Lanczos approach used here, which gives an 
important hopping contribution to the band gap. 

It is, nevertheless, true that hopping is suppressed in 
for large U, as one would expect for strongly correlated 
systems. The band gap is, however, an energy differ- 
ence between ground-state energies, where hopping has 
been suppressed differently. Thus hopping is suppressed 
very strongly in the TVM-particle state and less strongly 
in the (NM± l)-particle states, so that in the difference 
(Eq. (Q)) it appears as if the hopping had been enhanced. 

The reduced spins along the path of a moving extra 
occupancy plays two different roles. One role is that 
when the extra occupancy moves along different paths 
to a given site, the resulting states depend on this path. 
This is in contrast to the one-electron problem, and it 
has appreciable effects on the Lanczos matrix elements 
a n and 6 Tl +i. This is seen in the greater tendency to 
put weights into the diagonal elements a n in the one- 
particle case, and a reduction of the one-particle band 
width for non-bipartite lattices. For the non-bipartite 
clusters considered here, the hopping reduction of the 
band gap is therefore larger than the single-particle band 
width, even for N = 1, contrary to the intuitive feeling 
that the flipped spins should reduce the hopping con- 
tribution (see Table ||). The second effect is that the 
reduction of spins costs energy. This does not enter for 
U — > oo, since the energy cost (~ t 2 /U) for flipping the 
spin of one electron then goes to zero. Even for finite val- 
ues of U this does not enter in the first element at, but 
gradually becomes more important for the higher coeffi- 
cients. Normally, however, the spin flip energy is small 
compared with the band width W. Since the nondiago- 
nal elements b n +t are of the order \/NW or larger, these 
elements should tend to dominate the spin-flip energy, in 
particular in the large degeneracy limit. 
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